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1. Introduction 

Diffractive p meson photoproduction^ can be described within the dipole model [1-4] ac- 
cording to the formula 

^mAx{s, t)|,=o = sY^JA\dz vl/j^(r, z)a{s, r)^>''^{r, zf (1.1) 

h,h 

for the imaginary part of the forward elastic scattering amplitude, where ^1\ {r, z) and 

^{''j-(r, z) are the light-cone wavefunctions of the photon and vector meson, and a{s,r) is 
the dipole cross section. The light-cone wavefunctions represent the probability amplitudes 
for the photon or p meson to fluctuate into a qq pair of transverse size r in which the quark 
carries a fraction z of the meson's light-cone momentum. The sum is over quark/antiquark 
helicities {h and h) and X = L or T labels the polarization of the photon and meson. 
The photoproduction cross section is, after integrating over t, 

^xis) = ;^i^(3™-^a(^, 0)f (1 + /3f ) , (1.2) 

where /3a is the ratio of real to imaginary parts of the amplitude and B is the diffractive 
slope^. We shall assume that 
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^We use 'photoproduction' to refer to production using both real and virtual photons. 
■^We assume da/dt oc exp(_Bi). 



with A'^ = 0.55 GeV ^, which is in accord with the ZEUS data [5]. The HI data [6] prefer 
a somewhat larger value of B, but with a larger uncertainty. We compute /3\ according to 

/3a = tan (I^a) (1-4) 

with 

Qlnl^m^Al ., .. 

ax = oi /-I / ^ (1-5) 

aln(l/x) 

where Sm^A is given by Eq. (1.1) and x = (Q^ + 4?72^)/(Q^ + W'^) [nif is defined in the 
next section). The total cross section that is measured experimentally is given by 

a = aT + eaL (1.6) 

where e = 0.98. 

Our goal in this paper is to use the current HERA data on /o-meson photoproduction 
[5,6] to extract the meson's light-cone wavefunction. To do that we must specify the dipole 
cross section, which we do by assuming the FS2004 saturation model that was extracted 
from the HERA deep inelastic scattering data in Ref. [7]. 

2. Light-cone wavefunctions 

The photon's light-cone wavefunctions are [8-10]: 

*5(-' -) = \[§^K-mM^ - -)Q^ ' (2-1) 

and 

^lf{r, z) = ±/§ee; [ie^^'^ {z6,^-,^ - (1 - z)6,^-,^)dr + m;<5,±,,±] ^^ , (2.2) 

where 

e'^ = z{l - z)Q'^ + mj . (2.3) 

These wavefunctions are derived from perturbative QED and depend upon a phenomeno- 
logical light-quark mass, mf. We are compelled to take mj = 0.14 GeV, as determined by 
the fit in Ref. [7]. We take ej = l/V^, as appropriate for p meson production. 

The meson's light-cone wavefunctions can be written in terms of the scalar wavefunc- 
tions 4>l,t{'I', z): 

^U(-'-) = V^V-/^ M,z(l-z) ["^^ - ^)^P +"^/ - V?]'/'^(-'-) (2-4) 

where V^ ^ }:dr + d^ and 

^5(^' ^) = ^V^^(r^[^"^''''(^'^^±-^T - (1 - ^K^,-h±)dr + mf6^,^-^^]Mr, z). (2.5) 



These wavefunctions are subject to two important constraints. The first is the normaU- 
sation condition, which embodies our assumption that the p meson consists solely of a qq 
pair: 

j; /d2rdz|vl'j^^(r,z)|2 = l. (2.6) 

The second constraint arises from the measured value of fp, the meson decay constant for 
the longitudinally polarised meson, i.e. 

N r^ dz 
fpMp = -^ef / — -- [z{l-z)Ml + m)-VmL{r,z)l^^ . (2.7) 

" Jo ^{^ ^) 

The decay constant is deduced from the experimentally measured electronic decay width 
via the relation [10]: 






where T^^e+e- = 7.04 ± 0.06 keV [11]. 

3. Fitting the HERA data 

We must specify the form of the scalar wavefunctions (f)\. In Ref. [12], a 'Boosted Gaussian' 
(BG) wavefunction of the form 

2 d2 \ / ^2 d2 \ 
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4>f%,z) = AA, 4[z(l - z)r^2nRl exp ( ^j exp ( -^|^(Y^^ 1 (3-1) 

xexp^ ^2 

was used. This wavefunction is a simplified version of that proposed originally by Nemchik, 
Nikolaev, Predazzi and Zakharov [13]. 

In Ref. [12], it was assumed that bx = 1 and that Rl = Rt = R- The leptonic decay 
width constraint and the normalization conditions fix R and A/a (i.e. R^ = 12.3 GeV~^), 
leaving no free parameters. Predictions can then be made for the p-meson photoproduction 
cross section. This procedure leads to reasonable agreement with the old HERA data on 
the light vector mesons [12] and also for the heavier J/^ [14]. However, it is not able 
to accommodate the most recent HERA data on p production. Comparison with the 
HERA data leads to a x'^/data point of 234/75. For comparison with our later results, the 
longitudinal and transverse BG light-cone wavefunctions are shown in Figure 1. 

The poor agreement with data is considerably improved by allowing bx and Rx to vary 
freely.^ Specifically, we fit to the photoproduction total cross section,^ the longitudinal- 
to-transverse cross section ratio and to the electronic decay constant datum using Minuit 

^Allowing only Rx to vary does not much improve the fit. 

*We rescale the HI and ZEUS data by 0.95, which is consistent with the experimental uncertainty in 
the overall normalisation of the data [5,6]. 





Figure 1: The modulus squared of the longitudinal (left) and transverse (right) 'Boosted Gaussian' 
light-cone wavefunctions from Ref. [12]. 



[15]. To ensure that we fit in the diffractive region, we exclude those data points with 
T^ < 60 GeV. We fit to the most recent data from HI [6] and ZEUS [5]. In addition, we 
include the earlier HERA data [16-19] when they are in a kinematic region not covered 
by the latest data. This selected data set comprises 39 data points from ZEUS (32 points 
for the total cross section and 7 points for the longitudinal-to-transverse cross section 
ratio) and 36 data points from HI (27 points for the total cross section and 9 points for 
the longitudinal-to-transverse cross section ratio). To these, we add one data point for the 
decay constant of the longitudinally polarised p meson, giving a total of 76 data points. The 
result is a x^/degree of freedom equal to 82/72 and the corresponding best-fit parameters 
are listed in Table 1. The results of this fit can be seen as the dotted lines in Figures 2 and 
3, where the poor agreement at Q^ = is to be noted. 

Boosted Gaussian fit 



Free parameter 


Fitted value 


Rl 


27.33 GeV~2 


Hrp 


30.87 GeV-2 


hi 


0.5545 


br 


0.6792 



Table 1: The parameters corresponding to the BG fit of Eq. (3.1). The x^/degree of freedom is 
equal to 82/72. 



We can further improve the quality of fit (especially at Q^ = 0) by allowing for addi- 
tional end-point enhancement in the meson wavefunctions, i.e. using a scalar wavefunction 
of the form 

4>x{r, z) = 0f G(r, z) X [1 + cxe + dxt] (3.2) 



where ^ = 2z — 1. We find a preference for enhancement only in the wavefunction for 
transversely polarized mesons, i.e. the data prefer cl = cIl = 0, which is not surprising 
since Ui vanishes at Q^ = 0. The new fit lowers the x^/degree of freedom to 68/70. 
The resulting fitted parameters are listed in Table 2 and the corresponding cross section 
predictions appear as the solid lines in Figures 2 and 3. 

Best fit 



Free parameter 


Fitted value 


Rl 


26.76 GeV~2 


Rl 


27.52 GeV-2 


bL 


0.5665 


br 


0.7468 


Ct 


0.3317 


cIt 


1.310 



Table 2: The parameters corresponding to the BG fit with additional end-point enhancement in 
the transverse wavefunction, i.e. Eq. (3.2). The x^/degree of freedom is equal to 68/70. 



In Figure 4, we show the wavefunctions corresponding to the improved fit. Compared to 
Figure 1, the end-point enhancement is quite distinctive. Figure 5 shows the wavefunctions 
for both fits (the solid and dashed curves) at r = 0, plotted as a function of z. Note that 
they are almost indistinguishable from each other in the longitudinal case. For comparison, 
also shown on this figure is the result of the original BG parameterization [12]: it gives the 
dotted curves. 

The extent to which the data require any additional end-point enhancement should 
be set in context since the Q^ — t- limit suffers from the greatest theoretical uncertainty. 
For example, we characterize non-perturbative effects in the photon wavefunction through 
a single parameter (the quark mass) and it is unclear how a more sophisticated treatment 
would affect our conclusions. Moreover, it is also possible to improve the quality of fit to 
the data without appealing to Eq. (3.2) by increasing the value of the diffractive slope, B, 
at Q^ = 0. Specifically, we acheived a x^/degree of freedom of 67/72 after increasing B by 
15% at Q^ = relative to the value determined by Eq. (1.3). 

We note that both our fits indicate that the data prefer a transverse wavefunction with 
enhanced contributions at z — )• 0, 1. This conclusion is valid regardless of the uncertainties 
on the forward slope parameter. 

Finally, we should remark that it is possible to lower the x^ of the fit still further if 
we are prepared to rescale the data downwards by more than 5%. For example, rescaling 
down by 12% leads to a x^/degree of freedom equal to 62/70. The dashed curve in Figure 
6 confirms that the corresponding wavefunction is not very different from the one obtained 
by rescaling the data down by 5%. 



4. The leading-twist distribution amplitude 

We can use the longitudinal wavefunction from our fits to extract the corresponding leading 
twist-2 distribution amplitude, 93(2, /x) [20]: 

y,{z,fi)r^ jd^ke{\k\<fi)MKz) (4.1) 

where 

MKz)^ jd^re-'''-''Mr,z) , (4.2) 

which implies 

^{z,fj.) ~ / dr fiJi{nr)(j)L{r, z) . (4.3) 

Substituting for the scalar wavefunction from the previous section gives 

^{z, ^) ~ (1 - e-^'/^(^)') e-^'"//^(")' [z{l - z)]'^ , (4.4) 

where A(z)^ = 8[z(l — z)]^^/R'j^. We note that our DA is very slowly varying with ^ 
for fi > 1 GeV. This means that our parameterization neglects the perturbatively known 
//-dependence of the DA and can thus be viewed as a parameterization of the DA at some 
not too large a value of /x. This is reasonable given the limited Q^ range of the HERA data 
to which we fit (i.e. \/Q^ < 6 GeV). 

To compare with existing theoretical predictions for the DA, we compute moments of 
our DA, i.e. 



{C),= / dze>(z,/.). (4.5) 

Jo 

The n = moment is fixed by the decay constant constraint and is not a prediction^, we 
therefore follow convention and normalize the DA according to 

[ dz^{z,f,) = l. (4.6) 

Jo 

Our results are compared with the existing predictions in Table 3 and are in very good 
agreement with the expectations based on QCD sum rules and from the lattice. Also shown 
for comparison is the prediction based upon the old BG wavefunction used in Ref. [12], 
which does not fit the HERA data. The predictions in that case are rather similar to those 
of Ref. [21] using the light-front quark model. 

Finally, in Figure 7 we compare our DA with that predicted by Ball and Braun [22]. 
The agreement is reasonable given that in Ref. [22], the expansion in Gegenbauer polyno- 
mials is truncated at low order, which is presumably responsible for the local minimum 
at z = 1/2. Certainly the two distributions are broader than the asymptotic prediction 
~ 6z(l — z). 



^It is equivalent to Eq. (2.7) with the higher-twist nij — V^ terms set to zero. 



Moments of the leading twist DA at the scale jj, 



Reference 


Approach 


Scale /i 


{e). 


r)M 


r). 


{e). 


{e'). 


(This paper) 


Best fit 


~1 GeV 


0.227 


0.105 


0.062 


0.041 


0.029 


(This paper) 


BGfit 


~1 GeV 


0.229 


0.107 


0.063 


0.042 


0.030 


(This paper) 


Old EG prediction 


~1 GeV 


0.182 


0.072 


0.037 


0.022 


0.014 


[21] 


LFQM 


1 GeV 


0.19[0.21] 


0.08[0.09] 


0.04[0.05] 






[23] 


GcnSR 


1 GeV 


0.227(7) 


0.095(5) 


0.051(4) 


0.030(2) 


0.020(5) 


[24] 


SR 


1 GeV 


0.26 


0.15 








[22] 


SR 


1 GeV 


0.26(4) 










[25] 


SR 


1 GeV 


0.254 










[26] 


SR 


1 GeV 


n oo 1 0.03 
u.zoinQQ2 


1 1 +003 








[27] 


Lattice 


2 GeV 


0.24(4) 












6z{l - z) 


oo 


0.2 


0.086 


0.048 


0.030 


0.021 



Table 3: Our extracted values for (^");i, compared to predictions based on the light-front quark 
model (LFQM), QCD sum rules (SR), Generalised QCD Sum Rules (GenSR) or lattice QCD. Two 
predictions are given for each moment in the LFQM approach; one corresponding to a harmonic 
oscillator potential and the other (in square brackets) a linear potential. 

5. Conclusions 

The dipole model of diffractive photoproduction has been used successfully to describe a 
large body of data [14]. In this paper we have used it, together with accurate data on 
p- meson photoproduction collected at the HERA collider [5,6], in order to extract the p 
meson's light-cone wavefunction. The data require qualitatively different behaviour for the 
two meson polarizations. In particular, there is evidence for an enhancement of the end- 
point contributions to the wavefunction for transversely polarized mesons. We extracted 
the leading-twist /9-meson distribution amplitude for longitudinal polarization and found 
it to agree well with predictions based on QCD sum rules and from the lattice. 
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Figure 2: Comparison to the ZEUS data. The ctl/ctt data are at VF = 90 GeV. 
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Figure 3: Comparison to tlie HI data. Tlic olJut data are at W^ = 75 GeV. 
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Figure 4: The longitudinal (left) and transverse (right) light-cone wavefunctions squared corre- 
sponding to the BG fit with additional end-point enhancement in the transverse wavefunction. 
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Figure 5: The longitudinal (left) and transverse (right) light-cone wavefunctions squared corre- 
sponding to the BG fits with (solid) and without (dashed) additional end-point enhancement in the 
transverse wavefunction. The BG parameterization of Ref. [12] is also shown as the dotted curve. 
All curves evaluated at r = 0. 
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Figure 6: The longitudinal (left) and transverse (right) light-cone wavefunctions squared corre- 
sponding to the BG fits with additional end-point enhancement in the transverse wavefunction. The 
solid and dashed curves are extracted after rescaling the data down by 5% and 12% respectively. 
Both curves evaluated at r = 0. 
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Figure 7: The extracted DA at /^ = 1 GeV (sohd) compared to the DA at 1 GeV of Ref. [22] 
(dotted) and the asymptotic DA (dashed). 
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